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Abstract 

Recent analyses of the WMAP 5- year data constrain possible non-adiabatic contributions 
to the initial conditions of CMB anisotropies. Depending upon the early dynamics of the 
plasma, the amplitude of the entropic modes can experience a different suppression by the 
time of photon decoupling. Explicit examples of the latter observation are presented both 
analytically and numerically when the post-inflationary dynamics is dominated by a stiff 
contribution. 
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The initial conditions of the Boltzmann hierarchy can be usefully classified into adiabatic 
and non-adiabatic (i.e. entropic). The thermodynamic origin of this classification resides in 
the observation that the relative fluctuations of the specific entropy «Sy do not necessarily 
vanish (as in the case of the adiabatic mode). The relative fluctuations of the specific entropy 
<5j can be written, in gauge-invariant terms, as 

$, = ^U-3(Ci-<i), Ci = -* + -^T, (1) 

Sij Wi "T -L 

where W[ is the barotropic index of the i-th species; in Eq. (1) \1> represents the gauge- 
invariant Bardeen potential and 5i is the gauge-invariant density contrast. In the ACDM 
model (where A stands for the dark-energy component and CDM stands for the cold dark 
matter contribution) the indices i and j of Eq. (1) run over the four species of the plasma 
so that, in general, the initial conditions will contemplate one adiabatic mode and four non- 
adiabatic modes (see, for instance, [1, 2]). The initial conditions of the Boltzmann hierarchy 
can be set either by choosing only the adiabatic mode, or by selecting a combination of 
the adiabatic mode with one (or more) non-adiabatic modes. The obtained angular power 
spectra (both for temperature and polarization) can then be compared with the experimental 
data and interesting bounds can be set on the various combinations of the initial conditions 2 
[3, 4] (see also [1, 2]). 

Are there simple dynamical recipes able to suppress the entropic contributions? This 
is the basic question addressed in this paper. In the current framework, after inflation, 
the plasma was suddenly dominated by radiation. Absent the latter assumption, the non- 
adiabatic contribution to the pre-decoupling fluctuations of the spatial curvature will have a 
different relation to the entropic modes originally present right after inflation. While it is not 
mandatory to postulate different dynamical evolutions, it is useful to be aware of different 
possibilities which can help more dedicated scrutiny of the observational data. 

Prior to the radiation epoch, the plasma might have been expanding at a slower rate. This 
perspective was also invoked by Zeldovich who suggested that, prior to radiation dominance, 
the Universe was indeed quite stifff and characterized by a sound speed even comparable with 
the speed of light [5]. Post-inflationary phases stiffer than radiation can even lead to relic 
gravitons whose spectral energy density increases as a function of the comoving frequency 
[6]. If the inflaton field is identified with the quintessence field a stiff post-inflationary phase 
arises naturally [7] and this is what happens in the context of the so-called quintessential 
inflationary models [8] as well as in related contexts [9, 10]. 

Consider, for sake of simplicity, a post-inflationary plasma characterized by three distinct 
components so that the total energy density and the various pressures can be written as: 

Pt = Pm + Pr + PS, Pm = 0, p r = y p S = Wp S , (2) 

2 According to the present data, the adiabatic mode will have to be dominant in comparison with the 
remaining one (or more) entropic contributions. In the opposite case the (observed) anti-correlation peak in 
the temperature/polarization would not be correctly reproduced. 
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where p r and p m denote, respectively, the radiation and the matter energy densities while p r 
and p m indicate the corresponding pressures. In Eq. (2) ps and ps are the energy density 
and pressure of a supplementary component whose generic barotropic index will simply be 
denoted by w. From the pertinent Friedmann-Lemaitre equations 3 

n 2 = ^^p t , H 2 -H' = 4nGa 2 (p t + Pt ), H = -, (3) 

it can be easily argued that, for w > 1/3, ps dominates (at early times) in comparison with 
p m and p r . In Eqs. (3), ps — a~ 3 ^ w+1 ^ while p r ~ a~ 4 ; ergo ps will decrease faster than p r 
and this is the reason why the backreaction of the radiation can be important [6, 7, 8]. 

For the present purposes, it is practical to separate the adiabatic and the entropic fluc- 
tuations composing the gauge-invariant perturbation of the total pressure i.e. 



Spt = c 2 st 5p t + <5p nad , (4) 




where <kij are the fluctuations in the specific entropy already introduced in Eq. (1). In Eq. 
(5) the subscripts in the round brackets remind that the variation must be taken, respectively, 
for <5<Tij = and for Sp t = 0. Since, according to Eq. (2), the plasma is composed by three 
species there will be, in general, three entropic contributions. Equation (5) allows, indeed, 
for a more explicit form of Sp nSu d- 

SPnad = 2 APj(4 - C sj)^ij, 4 = ^, (6) 

0/ ~ l Pt ij Pi 

where the summation indices run over the three (or more) species of the fluid mixture. 
Equation (6) directly follows from Eqs. (4) and (5) by considering a generic pair of fluids 
and by summing over all the components. Using Eq. (2) into Eq. (6), 5p na d can be explicitly 
obtained: ^ 

^Pnad = ~ [PmPr^mr + 3wf/ m p' s S mS - (3«7 - l)f/ s p' T S Sl ], (7) 

where <S mr , S m $ and S$ r are, respectively, the three independent entropic fluctuations which 
can arise in the problem. To pass correctly from Eq. (6) to (7) it should be borne in mind 
that <Sij = — «Sji- In a democratic perspective all the entropic contributions in Eq. (7) can be 
present and with comparable amplitude. In the complementary situation one of the terms 
(e. g. <S mr ) is much larger than the remaining two. 

The fate of the non-adiabatic contributions given in Eq. (7) can be determined from 
the gauge-invariant evolution equations of the curvature and metric inhomogeneities. The 

3 A conformally flat metric g^ v = a 2 (r)?7^ will be assumed throughtout. The prime denotes a derivation 
with respect to the conformal time coordinate r. 
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gauge-invariant form of the Hamiltonian and the momentum constraints is 

V 2 ^ = AuGa 2 Pt e u V 2 (ft$ + *') = -AuGa\p t + p t )9 t , (8) 
(pt + Pt)0t = £(Pi + Pi)0i, (9) 

i 

where 9 t is the three-divergence of the (total) velocity field (defined in (9) as the sum of 
the velocities of the individual fluids) and e t is the gauge-invariant density contrast which 
corresponds to the total density contrasts in the comoving orthogonal gauge [13, 14, 15]. In 
terms of e t the momentum constraint (i.e. first relation of Eq. (8)) takes a form which is 
reminiscent of the (non-relativistic) Poisson equation. It turns out that e t is proportional to 
the difference of other two useful gauge-invariant quantities: 

Pt 

where 

c . E 4 Cl ,_*_^, "ffi+rv (id 

i Pt Pt A7vGa 2 (p t + p t ) 

The gauge-invariant variable ( can be interpreted either as the curvature perturbation in the 
uniform density gauge or as the density contrast in the uniform curvature gauge [13, 14, 15] 
(see also [16, 17]). The variable 1Z represent the (gauge-invariant) curvature perturbations 
which effectively correspond to the fluctuations of the spatial curvature on comoving orthog- 
onal hypersurfaces. The evolution of £ can be easily obtained from the equation for the total 
density fluctuation derived from the perturbation of the covariant conservation of the (total) 
energy-momentum tensor, i.e. 

5p[ - 3*'(pt + p t ) + (pt + fh)0t + 3n(5 Pt + S Pt ) = 0. (12) 

From Eq. (11) it can be easily deduced that 5p t = 3(p t + Pt)(C + ^)- Using the latter relation 
inside Eq. (12) the evolution of ( is simply 

C= ^Pnad-|, (13) 

Pt + Pt 3 

where Eqs. (3) and (4) have been used. Recalling Eqs. (10)— (11) the evolution equation for 
1Z can be directly obtained and it is 

n ' = X- ^ nad + 19 r v I ^ V2($ " " a rv° ! \ ^ 

p t + p t 127rGa 2 (p t + p t ) 47rGV(p t + p t ) 

The dynamical content of Eqs. (14) and Eq. (13) is clearly the same, i.e. the two equations 
coincide exactly in the long wavelength limit (i.e. kr <^ 1). Since p s ^ in Eq. (2), from 
Eqs. (7) and (14) the evolution for the curvature perturbations becomes then, to leading 
order in kr, 

W k = Tk(t,wJ s Jm), = -Qi(T,w,f s ,fM)^k ~ 02Cr,u>,/s,/ M )ftfc- (15) 
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Figure 1: The evolution of the curvature perturbations as a function of the scale factor and 
for various fixed values of w (plot at the left); the evolution of curvature perturbations for 
fixed a ^> 1 and as a function of the sound speed of the stiff component (plot at the right). 



In Eq. (15) the function J r n( T , w , fs, /m) is given by 

H{(w + l)[4(3w - lj/sW"" 1 - 9wf m f s S mS a 3w ] - 4/Aa 6 "- 1 } 
[4a 3 — 1 + 3a 3w f m + 3(w + l)/ s ] 2 

while the remaining two functions are: 



(16) 



r , f f s ^ [6^- 1 + (3 W + 5)/s + 5/ m a 3 -] 

ft(r, w, f s , fu) = H 2a[a ^ 1 + fs + fma3w] , (17) 

rf f f ^ ay [3/ m a 3 " , + 4a 3 - 1 + 3( W + l)/ s ] ns , 
&(r, /s, /m) = ft 2a[a3M> - 1 + /s + /ma3m] . (18) 

having introduced the following rescalings: 

a = — } f m = , / s = / x , a* = a(r*). (19) 

d* Pr(«*J Pr(«*J 

In what follows we will always take /s = 1 and set initial conditions for the numerical 
integration for ot\ = di/a* < 1. It is worth stressing that the condition /g = 1 implies 
that ps( a *) = Pr(o*). The value of / m < 1 specifies the fraction of non-relativistic matter 
eventually present for r*. The results of the numerical integration are reported in Figs. 1 
and 2 for some illustrative sets of parameters. With the full line (plot at the left in Fig. 
l)we have the result for the case S r s = <S m s = 0. In this case it is well known that the non- 
adiabatic solution to |72.jj.(a)| goes as <S mr /3 (units S mr = 1 will be adopted throughout). This 
is exactly the result of the full line reported in Fig. 1 (plot at the left) where, asyptotically 
for a ^> 1, log \TZk\ — —0.477. In the limit fs — > 0, the equation for IZk (see Eq. (15)) can 
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Figure 2: The evolutions of (plot at the left) and oiTZk (plot at the right) are illustrated. 
In the right plot the cases S m s = <S mr and iS m s <C S mr are compared. 

be written as 

(4 + 3f m a) 3 

where IZ^k) parametrizes the adiabatic solution which has been added for completeness but 
which will be left untouched by the present considerations. Equation (20) gives the value of 
the curvature perturbations induced by the non-adiabatic CDM-radiation mode when, right 
after inflation the Universe is dominated by a radiative equation of state. Recalling the form 
of (ordinary) Sachs- Wolfe contribution to the temperature anisotropies we will than have 
that, in the sudden decoupling approximation, A^ w ^ ~ —JZ*(k)/5 + 2<S*(&)/5. If the fluids 
do not exchange energy and momentum, as customarily assumed in the simplest situation 
and as verified in our case, then, S'^ = in the long wavelength limit. Indeed, it can be 
easily shown that 8[ = 3(wi + 1)^' — (wi + 1)9\ where 9\ = V • v\. Equation (1) then implies 
S[j = [6i — 0j), i.e. 5;j = up to corrections 0(k 2 r 2 ). 

When ps ^ 0, in the generic situation all the iSy are non-vanishing. Absent any specific 
knowledge of the initial conditions, the various entropic fluctuations can be expected to be 
comparable, i.e. S T s(k) ~ <Sms(k) ~ 5 m R. If we take, for instance, / m = 10 -2 (dashed line 
in the left plot) 1lk( a ) gets first to 1/3 (as implied by Eq. (20) in the absence of ps ) and 
then decreases by reaching, subsequently, a constant value which is between two and three 
orders of magnitude smaller than the putative asymptotic value which characterizes the case 
<S r s = 5 m s = (i.e. 1/3). As f m decreases the intermediate plateau get larger, and the 
asymptotic value gets progressively reduced. In the left plot of Fig. 1 it has been assumed 
that w — 1. In the plot at the right TZk{a) is reported for a fixed value of a (i.e., more 
specifically, «m = 10 9 ) but as a function of w. As w — > 1 the suppression can even be, 
depending on the parameters O(10~ 5 ). Of course the specific figure depends upon the other 
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parameters. At the same time it is clear that the amount of suppression depends upon the 
degree of stiffness, i.e. upon \w — 1/3|. 

The occurence that, for a while, IZk — (9(1/3) in Figs. 1 and 2 just means that, depending 
upon f m the terms containing /g can be neglected for intermediate values of a. The decrease 
in !Zk{oi) can be also explained. The analytical estimate can be separated into two steps, 
i.e. between «; and ax > 1 and between ax and «m ^> 1. In the first transition the terms 
proportional to f m can be neglected. The approximate result will then be: 

K k {a x ) ^ n*(k) + (^) - - + O(a^). (21) 

V3/ s /4a x +3(w + l)/ s V3/ s / 

For the second transition the relevant term will be the one proportional to S mr since the 
other terms can be neglected for a > 1 and for w > 1/3. By fixing the integration constant 
from Eq. (21) the asymptotic result (valid in the limit a ^> 1) 

ft fc (a M ) ~ 7^(fc) + — -——— ~ I — — + — . (22) 

3/s 3/ m a M + 4 V3/s 3 / 9a M /m 

If we do not assume a large hierarchy between Ss r and <S mr the first term in the second 
equality of Eq. (22) will approximately vanish and the second term will lead to a decrease 
of TZk{a). This is the effect observed in Fig. 1. If the values of the entropic fluctuations are 
drastically different (e.g. <S mr dominates against the others) then the results of Eq. (20) will 
be approximately true, at least asymptotically. This aspect is illustrated in Fig. 2 (plot at 
the right) for few cases where S mr is ten times larger than the other entropic contributions. 
Always in Fig. 2 (plot at the left) the evolution of ^ fc is reported. The full line in Fig. 2 
(plot at the left) corresponds to the case ps = where, following the same considerations of 
Eq. (20), ^k(oi) — Smr/5, i.e. log ^(a)! ~ —0.698 in units S mr = 1. The same patterns of 
suppressions discussed in the case of IZk also arise, as expected for Note that, finally, 
once TZk and are known to a given order in kr, the constraints of Eqs. (8)-(9) and (10) 
can be used to derive e t , C and t . 

The examples presented in this paper suggest that the entropic fluctuations can be dy- 
namically suppressed if, after inflation, there is a stiff contribution to the primeval plasma. 
The rationale for the obtained result depends both on a modification of the dynamics (at 
early time the stiff contribution dominates) and upon an interference effect between the vari- 
ous entropic contributions. In fine, the exercise presented here suggests that that the bounds 
on the non-adiabatic contribution obtained from the CMB data analysis do depend upon a 
number of specific assumptions on the early thermal history of the background geometry. 
In different terms, the entropy fluctuations used to set initial conditions of the Boltzmann 
hierarchy prior to equality might be already suppressed as a consequence of the preceding 
dynamical evolution. 
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